LibrerÃas
library(readxl)
library(tidyverse)
Load data
Data <- read_excel("Data.xlsx")
x <- "Año"
y <- c("Lat", "Long", "Altitud", "TMAX", "TMIN")
tabla_final <- data.frame(
"Variable"= character(),
"Trend" = numeric(),
"t" = numeric(),
"p" = numeric(),
"P95_max" = numeric(),
"P95_min" = numeric(),
"F" = numeric()
)
for (i in 1:5) {
tabla <- data.frame(
"Variable" = NA,
"Trend" = NA,
"t" = NA,
"p" = NA,
"P95_max" = NA,
"P95_min" = NA,
"F" = NA
)
tabla$Variable <- y[i]
model_g <- lm(formula(paste(y[i], paste(x, collapse = "+"), sep = " ~ ")), data = Data)
tabla$Trend <- model_g$coefficients[[2]]
tabla$t <- summary(model_g)$coefficients[2, 3]
tabla$p <- summary(model_g)$coefficients[2, 4]
tabla$X95_max <- confint(model_g, "Año", level = .95)[, 2]
tabla$X95_min <- confint(model_g, "Año", level = .95)[, 1]
tabla$F <- summary(model_g)$fstatistic[1]
tabla_final <- rbind(tabla_final, tabla)
}
DT::datatable(tabla_final,
class = "nowrap",
options = list(autoWidth = TRUE,
columnDefs = list(list(
width = "50%",
targets = 0))))
spp <- unique(Data$Especie)
compare.coeff <- function(b_g,se_g,b_i,se_i){
return((b_g-b_i)/sqrt(se_g^2+se_i^2))
}
tabla_ind <- data.frame(
"Spp" = character(),
"Variable"= character(),
"Trend" = numeric(),
"t" = numeric(),
"p" = numeric(),
"P95_max" = numeric(),
"P95_min" = numeric(),
"F" = numeric(),
"Dif" = numeric()
)
for (n in 1:length(spp)) {
ind <- filter(Data, Especie == spp[n])
if (nrow(ind) > 10) {
for (i in 1:5) {
tryCatch({
tabla <- data.frame(
"Spp" = NA,
"Variable" = NA,
"Trend" = NA,
"t" = NA,
"p" = NA,
"P95_max" = NA,
"P95_min" = NA,
"F" = NA,
"Dif" = NA
)
#General
model_g = lm(formula(paste(y[i], paste(
x, collapse = "+"
), sep = " ~ ")), data = Data)
tabla$Spp <- unique(ind[[1]])
tabla$Variable <- y[i]
model_i <-
lm(formula(paste(y[i], paste(
x, collapse = "+"
), sep = " ~ ")), data = ind)
tabla$Trend <- model_i$coefficients[[2]]
tabla$t <- summary(model_i)$coefficients[2, 3]
tabla$p <- summary(model_i)$coefficients[2, 4]
tabla$X95_max <- confint(model_i, "Año", level = .95)[, 2]
tabla$X95_min <- confint(model_i, "Año", level = .95)[, 1]
tabla$F <- summary(model_i)$fstatistic[1]
b_g <- summary(model_g)$coefficients[2,1]
se_g <- summary(model_g)$coefficients[2,2]
b_i <- summary(model_i)$coefficients[2,1]
se_i <- summary(model_i)$coefficients[2,2]
tabla$Dif <- 2*pnorm(-abs(compare.coeff(b_g,se_g,b_i,se_i)))
tabla_ind <- rbind(tabla_ind, tabla)
}, error = function(e) {
cat(
paste0("WARNING: Specie ", ind[1, 1], " variable (", y[i], ") has"),
conditionMessage(e),
"\n"
)
})
}
} else{
print(paste0("Data for ", ind[1, 1], " specie are insufficient"))
}
}
## [1] "Data for Aegialia arenaria specie are insufficient"
## WARNING: Specie Agolius abdominalis pecoudi variable (Altitud) has 0 (non-NA) cases
## [1] "Data for Agolius bonvouloiri bonvouloiri specie are insufficient"
## WARNING: Specie Agolius bonvouloiri cantabricus variable (Altitud) has 0 (non-NA) cases
## [1] "Data for Agrilinus sp. specie are insufficient"
## [1] "Data for Ahermodontus ambrosi specie are insufficient"
## WARNING: Specie Ahermodontus marini variable (Altitud) has 0 (non-NA) cases
## [1] "Data for Ammoecius brevis specie are insufficient"
## WARNING: Specie Ammoecius dentatus variable (Altitud) has subÃndice fuera de los lÃmites
## WARNING: Specie Brindalus porcicollis variable (Altitud) has subÃndice fuera de los lÃmites
## [1] "Data for Ceratophyus alloini specie are insufficient"
## [1] "Data for Ceratophyus holfmannseggi specie are insufficient"
## [1] "Data for Ceratophyus rossii specie are insufficient"
## [1] "Data for Cheironitis furcifer specie are insufficient"
## [1] "Data for Chilothorax brancoi specie are insufficient"
## WARNING: Specie Chilothorax conspurcatus variable (Altitud) has 0 (non-NA) cases
## [1] "Data for Chilothorax pictus specie are insufficient"
## [1] "Data for Erytus cognatus specie are insufficient"
## [1] "Data for Esymus sesquivittatus specie are insufficient"
## [1] "Data for Euoniticellus pallens specie are insufficient"
## [1] "Data for Heptaulacus algarbiensis specie are insufficient"
## [1] "Data for Heptaulacus brancoi specie are insufficient"
## [1] "Data for Jekelius bahilloi specie are insufficient"
## [1] "Data for Jekelius chersinus specie are insufficient"
## WARNING: Specie Jekelius hispanus variable (Altitud) has 0 (non-NA) cases
## WARNING: Specie Jekelius intermedius variable (Altitud) has subÃndice fuera de los lÃmites
## [1] "Data for Limarus zenkeri specie are insufficient"
## [1] "Data for Liothorax plagiatus specie are insufficient"
## [1] "Data for Melinopterus pubescens specie are insufficient"
## WARNING: Specie Neagolius heydeni variable (Altitud) has subÃndice fuera de los lÃmites
## [1] "Data for Neagolius montanus specie are insufficient"
## [1] "Data for Ochodaeus inermis specie are insufficient"
## [1] "Data for Odonteus armiger specie are insufficient"
## [1] "Data for Onthophagus albarracinus specie are insufficient"
## [1] "Data for Onthophagus merdarius specie are insufficient"
## [1] "Data for Onthophagus nuchicomis specie are insufficient"
## [1] "Data for Onthophagus semicomis specie are insufficient"
## WARNING: Specie Oromus alpinus variable (Altitud) has 0 (non-NA) cases
## [1] "Data for Oxyomus sylvestris specie are insufficient"
## [1] "Data for Parammoecius pyrenaeus specie are insufficient"
## [1] "Data for Parataenius simulator specie are insufficient"
## [1] "Data for Plagiogonus arenarius specie are insufficient"
## WARNING: Specie Planolinus fasciatus variable (Altitud) has subÃndice fuera de los lÃmites
## [1] "Data for Pleurophorus mediterranicus specie are insufficient"
## [1] "Data for Psammodius asper specie are insufficient"
## [1] "Data for Psammodius basalis specie are insufficient"
## [1] "Data for Psammodius laevipennis specie are insufficient"
## [1] "Data for Psammodius plicicollis specie are insufficient"
## WARNING: Specie Pseudacrossus sharpi variable (Altitud) has 0 (non-NA) cases
## [1] "Data for Rhyssemus convexus specie are insufficient"
## [1] "Data for Rhyssemus germanus specie are insufficient"
## WARNING: Specie Rhyssemus parallelus variable (Altitud) has subÃndice fuera de los lÃmites
## [1] "Data for Rhyssemus sulcatus specie are insufficient"
## [1] "Data for Scarabaeus pius specie are insufficient"
## WARNING: Specie Silphotrupes escorialensis opaculus variable (Altitud) has 0 (non-NA) cases
## WARNING: Specie Silphotrupes orocantabricus variable (Altitud) has 0 (non-NA) cases
## [1] "Data for Thorectes ferreri specie are insufficient"
## [1] "Data for Trypocopris vernalis specie are insufficient"
tabla_ind_lat <- filter(tabla_ind, Variable == "Lat")
DT::datatable(tabla_ind_lat,
class = "nowrap",
options = list(autoWidth = TRUE,
columnDefs = list(list(
width = "50%",
targets = 0))))
for ( i in 1:length(spp)){
ind <- filter(Data, Especie == spp[i])
if (nrow(ind) > 10) {
print(ggplot() +
geom_smooth(data= Data, aes(x = Año, y = Lat),col = "black", fill = "black", method = "lm") +
geom_smooth(data= ind, aes(x = Año, y = Lat),col = "red", fill = "red", method = "lm")+
ggtitle(paste0(spp[i]))+
labs(x= "Year", y = "Latitude")+
theme_minimal())
}else{
print(paste0("Data for ", ind[1, 1], " specie are insufficient"))
}
}
## [1] "Data for Aegialia arenaria specie are insufficient"
## [1] "Data for Agolius bonvouloiri bonvouloiri specie are insufficient"
## [1] "Data for Agrilinus sp. specie are insufficient"
## [1] "Data for Ahermodontus ambrosi specie are insufficient"
## [1] "Data for Ammoecius brevis specie are insufficient"
## [1] "Data for Ceratophyus alloini specie are insufficient"
## [1] "Data for Ceratophyus holfmannseggi specie are insufficient"
## [1] "Data for Ceratophyus rossii specie are insufficient"
## [1] "Data for Cheironitis furcifer specie are insufficient"
## [1] "Data for Chilothorax brancoi specie are insufficient"
## [1] "Data for Chilothorax pictus specie are insufficient"
## [1] "Data for Erytus cognatus specie are insufficient"
## [1] "Data for Esymus sesquivittatus specie are insufficient"
## [1] "Data for Euoniticellus pallens specie are insufficient"
## [1] "Data for Heptaulacus algarbiensis specie are insufficient"
## [1] "Data for Heptaulacus brancoi specie are insufficient"
## [1] "Data for Jekelius bahilloi specie are insufficient"
## [1] "Data for Jekelius chersinus specie are insufficient"
## [1] "Data for Limarus zenkeri specie are insufficient"
## [1] "Data for Liothorax plagiatus specie are insufficient"
## [1] "Data for Melinopterus pubescens specie are insufficient"
## [1] "Data for Neagolius montanus specie are insufficient"
## [1] "Data for Ochodaeus inermis specie are insufficient"
## [1] "Data for Odonteus armiger specie are insufficient"
## [1] "Data for Onthophagus albarracinus specie are insufficient"
## [1] "Data for Onthophagus merdarius specie are insufficient"
## [1] "Data for Onthophagus nuchicomis specie are insufficient"
## [1] "Data for Onthophagus semicomis specie are insufficient"
## [1] "Data for Oxyomus sylvestris specie are insufficient"
## [1] "Data for Parammoecius pyrenaeus specie are insufficient"
## [1] "Data for Parataenius simulator specie are insufficient"
## [1] "Data for Plagiogonus arenarius specie are insufficient"
## [1] "Data for Pleurophorus mediterranicus specie are insufficient"
## [1] "Data for Psammodius asper specie are insufficient"
## [1] "Data for Psammodius basalis specie are insufficient"
## [1] "Data for Psammodius laevipennis specie are insufficient"
## [1] "Data for Psammodius plicicollis specie are insufficient"
## [1] "Data for Rhyssemus convexus specie are insufficient"
## [1] "Data for Rhyssemus germanus specie are insufficient"
## [1] "Data for Rhyssemus sulcatus specie are insufficient"
## [1] "Data for Scarabaeus pius specie are insufficient"
## [1] "Data for Thorectes ferreri specie are insufficient"
## [1] "Data for Trypocopris vernalis specie are insufficient"
## Longitude ### Table
tabla_ind_lon <- filter(tabla_ind, Variable == "Long")
DT::datatable(tabla_ind_lon,
class = "nowrap",
options = list(autoWidth = TRUE,
columnDefs = list(list(
width = "50%",
targets = 0))))
for ( i in 1:length(spp)){
ind <- filter(Data, Especie == spp[i])
if (nrow(ind) > 10) {
print(ggplot() +
geom_smooth(data= Data, aes(x = Año, y = Long),col = "black", fill = "black", method = "lm") +
geom_smooth(data= ind, aes(x = Año, y = Long),col = "red", fill = "red", method = "lm")+
ggtitle(paste0(spp[i]))+
labs(x= "Year", y = "Longitude")+
theme_minimal())
}else{
print(paste0("Data for ", ind[1, 1], " specie are insufficient"))
}
}
## [1] "Data for Aegialia arenaria specie are insufficient"
## [1] "Data for Agolius bonvouloiri bonvouloiri specie are insufficient"
## [1] "Data for Agrilinus sp. specie are insufficient"
## [1] "Data for Ahermodontus ambrosi specie are insufficient"
## [1] "Data for Ammoecius brevis specie are insufficient"
## [1] "Data for Ceratophyus alloini specie are insufficient"
## [1] "Data for Ceratophyus holfmannseggi specie are insufficient"
## [1] "Data for Ceratophyus rossii specie are insufficient"
## [1] "Data for Cheironitis furcifer specie are insufficient"
## [1] "Data for Chilothorax brancoi specie are insufficient"
## [1] "Data for Chilothorax pictus specie are insufficient"
## [1] "Data for Erytus cognatus specie are insufficient"
## [1] "Data for Esymus sesquivittatus specie are insufficient"
## [1] "Data for Euoniticellus pallens specie are insufficient"
## [1] "Data for Heptaulacus algarbiensis specie are insufficient"
## [1] "Data for Heptaulacus brancoi specie are insufficient"
## [1] "Data for Jekelius bahilloi specie are insufficient"
## [1] "Data for Jekelius chersinus specie are insufficient"
## [1] "Data for Limarus zenkeri specie are insufficient"
## [1] "Data for Liothorax plagiatus specie are insufficient"
## [1] "Data for Melinopterus pubescens specie are insufficient"
## [1] "Data for Neagolius montanus specie are insufficient"
## [1] "Data for Ochodaeus inermis specie are insufficient"
## [1] "Data for Odonteus armiger specie are insufficient"
## [1] "Data for Onthophagus albarracinus specie are insufficient"
## [1] "Data for Onthophagus merdarius specie are insufficient"
## [1] "Data for Onthophagus nuchicomis specie are insufficient"
## [1] "Data for Onthophagus semicomis specie are insufficient"
## [1] "Data for Oxyomus sylvestris specie are insufficient"
## [1] "Data for Parammoecius pyrenaeus specie are insufficient"
## [1] "Data for Parataenius simulator specie are insufficient"
## [1] "Data for Plagiogonus arenarius specie are insufficient"
## [1] "Data for Pleurophorus mediterranicus specie are insufficient"
## [1] "Data for Psammodius asper specie are insufficient"
## [1] "Data for Psammodius basalis specie are insufficient"
## [1] "Data for Psammodius laevipennis specie are insufficient"
## [1] "Data for Psammodius plicicollis specie are insufficient"
## [1] "Data for Rhyssemus convexus specie are insufficient"
## [1] "Data for Rhyssemus germanus specie are insufficient"
## [1] "Data for Rhyssemus sulcatus specie are insufficient"
## [1] "Data for Scarabaeus pius specie are insufficient"
## [1] "Data for Thorectes ferreri specie are insufficient"
## [1] "Data for Trypocopris vernalis specie are insufficient"
## Elevation ### Table
tabla_ind_elev <- filter(tabla_ind, Variable == "Altitud")
DT::datatable(tabla_ind_elev,
class = "nowrap",
options = list(autoWidth = TRUE,
columnDefs = list(list(
width = "50%",
targets = 0))))
for ( i in 1:length(spp)){
ind <- filter(Data, Especie == spp[i])
if (nrow(ind) > 10) {
print(ggplot() +
geom_smooth(data= Data, aes(x = Año, y = Altitud),col = "black", fill = "black", method = "lm") +
geom_smooth(data= ind, aes(x = Año, y = Altitud),col = "red", fill = "red", method = "lm")+
ggtitle(paste0(spp[i]))+
labs(x= "Year", y = "Elevation")+
theme_minimal())
}else{
print(paste0("Data for ", ind[1, 1], " specie are insufficient"))
}
}
## [1] "Data for Aegialia arenaria specie are insufficient"
## [1] "Data for Agolius bonvouloiri bonvouloiri specie are insufficient"
## [1] "Data for Agrilinus sp. specie are insufficient"
## [1] "Data for Ahermodontus ambrosi specie are insufficient"
## [1] "Data for Ammoecius brevis specie are insufficient"
## [1] "Data for Ceratophyus alloini specie are insufficient"
## [1] "Data for Ceratophyus holfmannseggi specie are insufficient"
## [1] "Data for Ceratophyus rossii specie are insufficient"
## [1] "Data for Cheironitis furcifer specie are insufficient"
## [1] "Data for Chilothorax brancoi specie are insufficient"
## [1] "Data for Chilothorax pictus specie are insufficient"
## [1] "Data for Erytus cognatus specie are insufficient"
## [1] "Data for Esymus sesquivittatus specie are insufficient"
## [1] "Data for Euoniticellus pallens specie are insufficient"
## [1] "Data for Heptaulacus algarbiensis specie are insufficient"
## [1] "Data for Heptaulacus brancoi specie are insufficient"
## [1] "Data for Jekelius bahilloi specie are insufficient"
## [1] "Data for Jekelius chersinus specie are insufficient"
## [1] "Data for Limarus zenkeri specie are insufficient"
## [1] "Data for Liothorax plagiatus specie are insufficient"
## [1] "Data for Melinopterus pubescens specie are insufficient"
## [1] "Data for Neagolius montanus specie are insufficient"
## [1] "Data for Ochodaeus inermis specie are insufficient"
## [1] "Data for Odonteus armiger specie are insufficient"
## [1] "Data for Onthophagus albarracinus specie are insufficient"
## [1] "Data for Onthophagus merdarius specie are insufficient"
## [1] "Data for Onthophagus nuchicomis specie are insufficient"
## [1] "Data for Onthophagus semicomis specie are insufficient"
## [1] "Data for Oxyomus sylvestris specie are insufficient"
## [1] "Data for Parammoecius pyrenaeus specie are insufficient"
## [1] "Data for Parataenius simulator specie are insufficient"
## [1] "Data for Plagiogonus arenarius specie are insufficient"
## [1] "Data for Pleurophorus mediterranicus specie are insufficient"
## [1] "Data for Psammodius asper specie are insufficient"
## [1] "Data for Psammodius basalis specie are insufficient"
## [1] "Data for Psammodius laevipennis specie are insufficient"
## [1] "Data for Psammodius plicicollis specie are insufficient"
## [1] "Data for Rhyssemus convexus specie are insufficient"
## [1] "Data for Rhyssemus germanus specie are insufficient"
## [1] "Data for Rhyssemus sulcatus specie are insufficient"
## [1] "Data for Scarabaeus pius specie are insufficient"
## [1] "Data for Thorectes ferreri specie are insufficient"
## [1] "Data for Trypocopris vernalis specie are insufficient"
## Maximum temperature ### Table
tabla_ind_tmax <- filter(tabla_ind, Variable == "TMAX")
DT::datatable(tabla_ind_tmax,
class = "nowrap",
options = list(autoWidth = TRUE,
columnDefs = list(list(
width = "50%",
targets = 0))))
for ( i in 1:length(spp)){
ind <- filter(Data, Especie == spp[i])
if (nrow(ind) > 10) {
print(ggplot() +
geom_smooth(data= Data, aes(x = Año, y = TMAX),col = "black", fill = "black", method = "lm") +
geom_smooth(data= ind, aes(x = Año, y = TMAX),col = "red", fill = "red", method = "lm")+
ggtitle(paste0(spp[i]))+
labs(x= "Year", y = "Maximum temperature")+
theme_minimal())
}else{
print(paste0("Data for ", ind[1, 1], " specie are insufficient"))
}
}
## [1] "Data for Aegialia arenaria specie are insufficient"
## [1] "Data for Agolius bonvouloiri bonvouloiri specie are insufficient"
## [1] "Data for Agrilinus sp. specie are insufficient"
## [1] "Data for Ahermodontus ambrosi specie are insufficient"
## [1] "Data for Ammoecius brevis specie are insufficient"
## [1] "Data for Ceratophyus alloini specie are insufficient"
## [1] "Data for Ceratophyus holfmannseggi specie are insufficient"
## [1] "Data for Ceratophyus rossii specie are insufficient"
## [1] "Data for Cheironitis furcifer specie are insufficient"
## [1] "Data for Chilothorax brancoi specie are insufficient"
## [1] "Data for Chilothorax pictus specie are insufficient"
## [1] "Data for Erytus cognatus specie are insufficient"
## [1] "Data for Esymus sesquivittatus specie are insufficient"
## [1] "Data for Euoniticellus pallens specie are insufficient"
## [1] "Data for Heptaulacus algarbiensis specie are insufficient"
## [1] "Data for Heptaulacus brancoi specie are insufficient"
## [1] "Data for Jekelius bahilloi specie are insufficient"
## [1] "Data for Jekelius chersinus specie are insufficient"
## [1] "Data for Limarus zenkeri specie are insufficient"
## [1] "Data for Liothorax plagiatus specie are insufficient"
## [1] "Data for Melinopterus pubescens specie are insufficient"
## [1] "Data for Neagolius montanus specie are insufficient"
## [1] "Data for Ochodaeus inermis specie are insufficient"
## [1] "Data for Odonteus armiger specie are insufficient"
## [1] "Data for Onthophagus albarracinus specie are insufficient"
## [1] "Data for Onthophagus merdarius specie are insufficient"
## [1] "Data for Onthophagus nuchicomis specie are insufficient"
## [1] "Data for Onthophagus semicomis specie are insufficient"
## [1] "Data for Oxyomus sylvestris specie are insufficient"
## [1] "Data for Parammoecius pyrenaeus specie are insufficient"
## [1] "Data for Parataenius simulator specie are insufficient"
## [1] "Data for Plagiogonus arenarius specie are insufficient"
## [1] "Data for Pleurophorus mediterranicus specie are insufficient"
## [1] "Data for Psammodius asper specie are insufficient"
## [1] "Data for Psammodius basalis specie are insufficient"
## [1] "Data for Psammodius laevipennis specie are insufficient"
## [1] "Data for Psammodius plicicollis specie are insufficient"
## [1] "Data for Rhyssemus convexus specie are insufficient"
## [1] "Data for Rhyssemus germanus specie are insufficient"
## [1] "Data for Rhyssemus sulcatus specie are insufficient"
## [1] "Data for Scarabaeus pius specie are insufficient"
## [1] "Data for Thorectes ferreri specie are insufficient"
## [1] "Data for Trypocopris vernalis specie are insufficient"
## Minimum temperature ### Table
tabla_ind_tmin <- filter(tabla_ind, Variable == "TMIN")
DT::datatable(tabla_ind_lat,
class = "nowrap",
options = list(autoWidth = TRUE,
columnDefs = list(list(
width = "50%",
targets = 0))))
for ( i in 1:length(spp)){
ind <- filter(Data, Especie == spp[i])
if (nrow(ind) > 10) {
print(ggplot() +
geom_smooth(data= Data, aes(x = Año, y = TMIN),col = "black", fill = "black", method = "lm") +
geom_smooth(data= ind, aes(x = Año, y = TMIN),col = "red", fill = "red", method = "lm")+
ggtitle(paste0(spp[i]))+
labs(x= "Year", y = "Minimum temperature")+
theme_minimal())
}else{
print(paste0("Data for ", ind[1, 1], " specie are insufficient"))
}
}
## [1] "Data for Aegialia arenaria specie are insufficient"
## [1] "Data for Agolius bonvouloiri bonvouloiri specie are insufficient"
## [1] "Data for Agrilinus sp. specie are insufficient"
## [1] "Data for Ahermodontus ambrosi specie are insufficient"
## [1] "Data for Ammoecius brevis specie are insufficient"
## [1] "Data for Ceratophyus alloini specie are insufficient"
## [1] "Data for Ceratophyus holfmannseggi specie are insufficient"
## [1] "Data for Ceratophyus rossii specie are insufficient"
## [1] "Data for Cheironitis furcifer specie are insufficient"
## [1] "Data for Chilothorax brancoi specie are insufficient"
## [1] "Data for Chilothorax pictus specie are insufficient"
## [1] "Data for Erytus cognatus specie are insufficient"
## [1] "Data for Esymus sesquivittatus specie are insufficient"
## [1] "Data for Euoniticellus pallens specie are insufficient"
## [1] "Data for Heptaulacus algarbiensis specie are insufficient"
## [1] "Data for Heptaulacus brancoi specie are insufficient"
## [1] "Data for Jekelius bahilloi specie are insufficient"
## [1] "Data for Jekelius chersinus specie are insufficient"
## [1] "Data for Limarus zenkeri specie are insufficient"
## [1] "Data for Liothorax plagiatus specie are insufficient"
## [1] "Data for Melinopterus pubescens specie are insufficient"
## [1] "Data for Neagolius montanus specie are insufficient"
## [1] "Data for Ochodaeus inermis specie are insufficient"
## [1] "Data for Odonteus armiger specie are insufficient"
## [1] "Data for Onthophagus albarracinus specie are insufficient"
## [1] "Data for Onthophagus merdarius specie are insufficient"
## [1] "Data for Onthophagus nuchicomis specie are insufficient"
## [1] "Data for Onthophagus semicomis specie are insufficient"
## [1] "Data for Oxyomus sylvestris specie are insufficient"
## [1] "Data for Parammoecius pyrenaeus specie are insufficient"
## [1] "Data for Parataenius simulator specie are insufficient"
## [1] "Data for Plagiogonus arenarius specie are insufficient"
## [1] "Data for Pleurophorus mediterranicus specie are insufficient"
## [1] "Data for Psammodius asper specie are insufficient"
## [1] "Data for Psammodius basalis specie are insufficient"
## [1] "Data for Psammodius laevipennis specie are insufficient"
## [1] "Data for Psammodius plicicollis specie are insufficient"
## [1] "Data for Rhyssemus convexus specie are insufficient"
## [1] "Data for Rhyssemus germanus specie are insufficient"
## [1] "Data for Rhyssemus sulcatus specie are insufficient"
## [1] "Data for Scarabaeus pius specie are insufficient"
## [1] "Data for Thorectes ferreri specie are insufficient"
## [1] "Data for Trypocopris vernalis specie are insufficient"